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ABSTRACT 

Context. The emission nebula N66 is the brightest H II region in the Small Magellanic Cloud (SMC), the stellar association NGC 346 
being located at its center. The youthfulness of the region NGC 346/N66 is well documented by studies of the gas and dust emission, 
and the detection in the optical of a rich sample of pre-main sequence (PMS) stars, and in the mid- and far-IR of young stellar objects 
(YSOs). However, a comprehensive study of this region has not been performed in the near-IR that would bridge the gap between 
previous surveys. 

Aims. We perform a photometric analysis on deep, seeing-limited near-IR VLT images of the region NGC 346/N66 and a nearby con- 
trol field of the SMC to locate the centers of active high- and intermediate-mass star formation by identifying near-IR bright objects 
as candidate stellar sources under formation. 

Methods. We use archival imaging data obtained with the high-resolution camera ISAAC at VLT of NGC 346/N66 to construct the 
near-IR color-magnitude (CMD) and color-color diagrams (C-CD) of all detected sources. We investigate the nature of all stellar 
populations in the observed CMDs, and we identify all stellar sources that show significant near-IR excess emission in the observed 
C-CD. We, thus, select the most likely young stellar sources. 

Results. Based on their near-IR colors, we select 263 candidate young stellar sources. This sample comprises a variety of objects, 
such as intermediate-mass PMS and Herbig Ae/Be stars and possible massive YSOs, providing original near-IR colors for each of 
them. The spatial distribution of the selected candidate sources indicates that they are located along the dusty filamentary structures 
of N66 seen in mid- and far-IR dust emission and agrees very well with that of previously detected candidate YSOs and PMS stars. 
Conclusions. Our study provides an original accurate set of near-IR colors for candidate young stellar sources. This provides signifi- 
cant information about the star formation process in NGC 346/N66, but does not establish the types of these objects, which requires the 
construction of complete spectral energy distributions for individual sources from multiwavelength data. This would be an important 
follow-up study to that presented here. 

Key words, stars: formation, pre-main-sequence - Magellanic Clouds - open clusters and associations: individual: NGC 346 - HII 
regions - ISM: individual objects: LHA 1 15-N66 

1. Introduction dHilditchet al.ll2005l) . It is embedded in the brightest H n re- 

,. , „ „ . gion of the SM C, which is referred to as LHA 1 15-N66 or N66 

Located in the stellar constellation il ucana, the Small Magellanic denize! [T956). With an Ha luminosity almost 60 times higher 

Cloud (SMC), is an excellent laboratory for investigating the ^ ^ star _ formi ion of orion ( K ennicutt 1984), N66 has 

star formation processes and the associated chemical evolu- a diameter of about ? , co ondi to ap p roxirriate ly 123 pc. 

tion of dwarf galaxies. Its present subsolar chemical abundance „, „ . XT ^^ o^/xt^ 

m nnn/i onn? c i \ ■ i- **u- 1 u The star-lorming region NGC 346/N66 comprises a van- 

(Z = 0.004; ~ 20% or solar) implies that this galaxy may have . „ , • — ,'■ ■ — —r\\ ~ nn A > r™ 

, . . . ., . . ,. ° £ ./ i . ety or young stellar populations Gouliermis et al. 2006). The 

characteristics similar to those in earlier times or the evolution J . . ^t-^A^a^ , — TT, 'T — 

„. „,,„. , „ , m_ stellar association NGC 346, located at the center or the so- 

or the universe. Since SMC is so close to our Galaxy, it is there- „ , . , „ , „ XT ,, , , , , - _ 

, „ . , , , „ j., , called nebular bar or N66, hosts the largest sample or O- 

rore an excellent laboratory among the large collection or dwarr . , „, fr-= nrrT^ i „, „ n 

, , , , \ . B . t °, . , f , , type stars in the entire SMC Massev et al. 1989; Walborn et al. 

irregulars and blue compact galaxies lor the study of resolved uVEwJ c; ^ , 1^™A ^ J"- — Z — j 1 ■ — : ^ 

w i j. - - rr, 2000; Evans et al. 2006). Studies based on deep imaging with the 

extragalactic stellar populations and star-forming regions. The ' " — ~ zrf , . . . „ f . . 

young stellar association NGC 346 (RA (J2000) = 00 h 59 m 18 s , Huhhle f *^S^ C that .*^f. V1 f imty °! the Wh °^ re " 

DEC (J2000) = -72° 10' 48") is a large star-forming cluster 8 10n °l^ C 346/N66 18 nc J.?» loW " and ln ™ dlate " 



in the SMC, located at a distance of about 60.6 kpc from us 



mass PMS stars, some of which exhibit recent star formation 
with an age $ 5 Myr, while others belong to an older underly- 
Send offprint requests to: D. A. G. e-mail: dgoulier@mpia.de ing_population of age < 10 Myr (see e.g.,|Hennekemperet al.| 

* Based on observations made with ESO Telescopes at the La Silla EQQS». Additional evidence of the youthfulness of this region 
Paranal Observatory under program ID 063.1-0329. comes from observations with the Spitzer Space Telescope 

** Table [TJ is available in its entirety only in electronic form at the and the discovery of 111 candidate massi ve young stellar ob - 
CDS via anonymous ftp to [f tp://cdsarc.u-strasbg T&1 or via jects (YSOs) with 2 < M/M < 17 dSimon et al] 12007). 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/515/A56 Nevertheless, although these objects should emit in near-IR 
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Fig. 1. Mosaic images from the combination of the northern and center jittered ISAAC frames of the area of NGC 346/N 66 in the 
(a) J-, (b) H-, and (c) K s -band. 



bands, no detailed study in such wavelengths exists in the lit- 
erature. 

Pioneering work on the dus t and gas conten t of 
NGC 346/N6 6 was performed by iRubio et all d2000l) and 
IContursi et alJ (120001) . who found a correlation between H2 
infrared emission and CO lines, characteristic of a photo- 
dissociated region (PDR). A PDR develops when the far-UV 
radiation of the bright OB stars reaches the surface of the 
parental molecular cloud. The degree of ionization decreases 
outwards, and a thin barrier develops that segregates the ionized 
from the atomic gas. While H2 is not fully ionized behind this 
front but partly dissociated, CO molecules located little deeper 
in th e cloud are more ea sily dissociated by absorbing UV pho- 
tons. |Rubio!etil] (|2003) and lContursi et all d2000l) inferred that 
star formation in NGC 346/N66 has taken place as a sequential 
process in the "bar" of N 66, which these authors define as 
the oblique bright emission region extending from southeast to 
northwest centered on NGC 346. This process results in several 
embedded sources, seen as IR-emission peaks in the 2.14 yum 
H 2 line and the ISOCAM LW2 band (5 - 8 urn). These peaks 
are alphabetically nu mbered from "A" to "I" ( Rubio et alJ l2000: 
IContursi et alj|20doh . the association NGC 346 itself coinciding 
with peak "C". 

Studies with the Advanced Camera for Surveys (ACS) on- 
board the Hubble Space Telescope (HST) have shown that the 
PMS stellar content of NGC 346/N66 covers a mass range in the 
subsolar regime. These studies sug gest that recent star formation 
occur red around 3-5 Myr ago dNota et all l2006t ISabbi et all 
l2007h there being an underlying older PMS population, which 
indicates that there were earlier star formation even ts that oc- 
curred about 10 Myr ago dHennekemper et alJl2008h . The PMS 
population is mainly centrally concentrated apart from the as- 
sociation NGC 346 in a number of subclusters with clustering 
proper ties quite similar to t hose of Milky Way star-forming re- 
gions (Schmei a et aL II2OO9I). Furthermore, the intermediate-age 
star clust er BS 90 dBica & Schmittlll995h with an age of about 
4.5 Gyr dRochau et al.ll2007h is also located in this part of the 
SMC, projected in front of N66. 

This photometric study of the star-forming region 
NGC 346/N66 focuses on identifying young stellar sources, 



mainly intermediate- and high-mass PMS stars of the region, 
stars that have not yet started their lives on the main sequence. 
The PMS phase in the evolution of stars with masses up to 
~ 6 M Q corresponds to the time between the gravitational core 
collapse, which forms the protostar (on the birthline), and the 
ignition of hydrogen in the formed star, placing it on the Zero- 
Age Main Sequence (ZAMS). During this evolutionary phase, 
the observed radiation from the star is affected significantly by 
circumstellar disks of dust and gas, formed by matter infalling 
during the collapse of the rotating core, and surface activity. 
Examples of these typical PMS stars, the T Tauri stars, thus, 
exhibit prominent optical emission lines, which are understood 
to stem from chromospheric heating, periodic fluctuations in 
light that indicate rotating star-spots, variability, and excess 
broadband flux in UV and IR, and are sometimes associated with 
molec ular outflows, winds, or accretion (e.g., Lada & Kvlafis 
119911) . Intermediate-mass (2 < M/M Q < 10) PMS stars ar e 
called Herbig AeBe (HAeBe) stars (e.g.. lPerez & Gradvl[T997h . 
Being more massive analogues of T Tauri stars, they are PMS 
A- and B-type stars that exhibit emission lines produced by 
both strong stellar winds and the cocoons of remnant gas from 
which they collapsed. They typically contain circumstellar 
disks and therefore have spectral energy distributions (SEDs) 
of young stellar objects (YSOs) of class II. They range in age 
between 0.5 and 5 Myr, similar to T Tauri stars. Because of the 
aforementioned characteristics, PMS stars are visible in near-IR 
wavebands. 

While extensive observations have been performe d of young 



stellar sources in the region NGC 346/N66 at optical (Nota et al 
2006) and mid- and far-IR (Bola tto et al 



2006 



Gouliermis et al 



2007: iSimon et al.l l2007) wavelengths, there has been no com- 
parative investigation in the near-IR. The present study aims to 
close this gap in the available spectral coverage for this region 
by acquiring near-IR data to characterize with greater accuracy 
its young stellar population. This will also aid our understand- 
ing of star formation in the low-metallicity environment of the 
SMC. In this paper, we present our ground-based near-IR pho- 
tometry derived from observations with VLT/ISAAC of the re- 
gion NGC 346/N66. We present the observational material used 
and its reduction in Sect. 2, and we discuss the photometric pro- 
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Fig. 2. Typical photometric uncertainties in the three /, H, and K s bands, derived by DAOPHOT for all detected sources. 



cess in Sect. 3. The various observed stellar populations and the 
corresponding stellar systems comprised in the observed field 
are discussed in terms of variations in the constructed color- 
magnitude diagrams in Sect. 4. We present the constructed near- 
IR color-color diagram of all detected sources and apply a se- 
lection criterion by identifying stars currently forming by means 
of their near-IR excess emission in terms of their positions in 
this diagram in Sect. 5. In the same section, the spatial distribu- 
tion and the nature of the selected young stellar sources is also 
discussed. Finally, conclusive remarks on this study are given in 
Sect. 6. 



2. Observations and data reduction 

Red giants, stars close the end of their life usually have high 
mass-loss rates, while stars at the earliest stages of their forma- 
tion are embedded into dense molecular clumps and cores. In 
both cases, the dust shells or disks around the objects of interest 
absorb almost all of the visible radiation, which is re-radiated 
at longer wavelengths. As a consequence, ab sorption decreases 
very r apidly with increasing wavelength (e. g jjovcel 1 992t TGlassI 
1999), the extinction coefficient at 2.2 yum being approximately 
10% of that at 500 nm. Here we are interested in young stellar 
sources associated with circumstellar dusty shells or disks, char- 
acterized by bright IR excess emission. 




2.1. Observations 

The near-IR images of NGC 346/N66 are obtained within the 
ESO Program ID 63.1-0329 (PI: M. Rubio) with the Infrared 
Spectrometer And Array Camera (ISAAC), mounted on the 
Nasmyth-Focus B of UT1 (Antu) at the Very Large Telescope 
(VLT). The images were taken between 24 and 26 Sep. 1999, 
under fair seeing conditions (FWHM between 0.6" and 1.0"). 
They were used p artly to detect YSOs in the region of N66 
dSimon et al. 2007), but have never been presented in their com- 
plete spatial and wavelength coverage. The observations were 
performed with the short wavelength arm of ISAAC, using the 
1024 x 1024 HgCdTe Hawaii Rockwell array. The pixel scale 
of the Hawaii detector is 0.1484"/pixel, providing a maximum 
field-of-view of 152 x 152 arcsec 2 . The images were obtained in 
the filters J (1.25 ^m), H (1.65 jum), and K, (2.16 fim). 

At wavelengths longer than ~ 2.2 /urn, the thermal back- 
ground is dominated by atmospheric and telescope emission, 
leading to a highly variable sky brightness in the infrared. 
Detector cosmetics and instabilities also illustrate the need for an 
accurate sky subtraction. This is achieved by mean of the tech- 
nique of jittering which is available for the short-wavelength 
and some long-wavelength imaging templates of ISAAC. A set 
of 10 to 100 frames are combined to form one final "jittered" 
frame. At the beginning and end of each night, twilight flats are 
taken for each filter. The following morning, dark images of the 
detector are taken for all detector integration times (DITs) used 
during the previous night. The area of NGC 346 was divided into 
two jittered frames, each consisting of 10 single frames. An ad- 
ditional jittered frame was taken for a control field of the galaxy, 
located to the south of the system. For each single frame, the to- 
tal integration time was 60s, taken in 5s (DIT) exposures of 12 
sub-integrations (NDIT). The field-of-view is about 2.5' x 2.5', 
correspondin g to about 44 x 44 p c 2 at the SMC (distance from 
us 60.6 kpc: iHilditch et alj 12005b . We retrieved the dataset of 
these observations, including the science frames, twilight flats, 
dark images, and standard stars from the ESO Science Archive 
Facilit>B 



2.2. Data reduction 

To reduce the data, we used ESO's ISAAC pipeline, which is 
based entirely on the libraries of the data reduction package 
ECLIPSE0. The package consists of the main eclipse-library and 
the additional instrument specific pipeline packages. Another al- 



1B 18 20 

apparent magnitude 



Fig. 3. Completeness for all stars detected in the three observed 
fields in all three J, H and K s filters. 



1 http : //www. eso . org/projects/dfs/papers/j itter99/ 

2 |http : / / archive . eso . org/ cms/ 1 

3 |http : //www . eso . org/eclipse] 
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gorithm of ECLIPSE is the jitter routine, which implements ef- 
ficient filtering and processing methods for infrared data reduc- 
tion. The detector pixels have intrinsically different sensitivities 
because of the quantum efficiency variations from pixel to pixel. 
There is also a readout gradient of the array, which itself is not 
homogeneously illuminated by the telescope. These effects, as 
well as artifacts caused by dust on the optical surface, are cor- 
rected within the data reduction process. 

The dark current of the ISAAC Hawaii array is low, so the 
detector bias, also called the "zero level offset", is the domi- 
nant feature of the dark frames. A master dark, constructed with 
ECLIPSE, was subtracted from the science frames. The effects 
of dust, the inhomogeneity of the illumination, and pixel to pixel 
variations were then removed using a master Hat-field frame. 
Since no flat-field screen exists at UT1, for ISAAC only twi- 
light flats are available. For the subtraction of both sky and bias, 
the darkest twilight flat frame was used. After this subtraction, 
the remaining 15 twilight flats were combined to produce a mas- 
ter flat-field with ECLIPSE, which also creates a bad pixel map 
output file. The jitter routine in ECLIPSE was then used to re- 
duce the science data for the frames observed with acceptable 
seeing, by filtering out low-frequency sky variations from the 
set of jittered images, a method called sky combination. Two 
sets of jittered images were combined for the whole area of 
NGC 346. The corresponding final frames are shown for each 
filter in Fig. Q] The association NGC 346, the numerous young 
compact star clusters located in the bright Hn-region N 66, and 
the intermediate-age star cluster BS 90 are easily identified in 
these images. The photometry, discussed in the following sec- 
tion, was performed on the single jittered frames for each filter. 
The derived catalogs were then matched to produce the complete 
photometric catalog of the observed regions. 



3. Photometry 

Our photometry was performed within the Image Reduction and 
Analysis Facil ity (IRAF) systerrQ, with the package DAOPHOI0 
(Stets onll987l) . The effect of both the quality of the imaging sys- 
tem and the seeing conditions on the obtained magnitudes of the 
stars is quantified by the point spread function (PSF), which rep- 
resents the image of a typical point source in the observations. 
DAOPHOT provides accurate PSF photometry of crowded stel- 
lar fields. The magnitudes of all stars detected in the observed 
jittered images are defined relative to the brightness of Vega. We 
first have to specify the input parameters of DAOPHOT to opti- 
mize the detection and photometry of stars, before using the rou- 
tine daofind to detect all stellar sources in the science frames. 
The routine phot is then used to perform aperture photometry of 
the detected stars and determine their instrumental magnitudes. 
The radius of the aperture is selected to be between 1 and 2 times 
that of the FWHM, depending on the effective seeing. For indi- 
vidual stars, the aperture radius can be up to 5 or more of that of 
the FWHM. 

The typical PSF of our images was constructed by interac- 
tively selecting approximately 20 isolated bright stars per frame 
using DAOPHOT pstselect. The PSF was modelled to be the 
sum of an analytic bivariate Gaussian function and empirical 
corrections from the best Gaussian of the true observed bright- 
ness values within the average profile of several stars in the im- 
age. This process and the computation of the PSF to be fitted 
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Fig. 4. Chart of all stars detected with good photometric ac- 
curacy in all three filters in all three fields observed with 
VLT/ISAAC. This map shows the relative positions of the ob- 
served fields: The northern field that includes the intermediate- 
age cluster BS 90, the central field, where the association 
NGC 346 and the main part of the H II region N66 are lo- 
cated, and the remote control-field, which covers the represen- 
tative stellar populations of the general SMC field in the area. 
North is up and East is to the left of the map. 



was performed by the routine psf. PSF photometry was per- 
formed with the routine allstar, which after classifying the 
stars in groups, compiles a catalog of the most likely candidate 
stars, based on their PSF fitting and the physical conditions, and 
subtracts them from the original image. The photometric pro- 
cess was repeated for the subtracted frames and the magnitudes 
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Table 1. Sample of the photometric catalog of stars found in this study in the region of NGC 346/N66 and a nearby control field in 
all three 7, H, and K s filters with VLT ISAAC imaging. This table is available in its entirety at the CDS. 



ID 


R.A. 


Decl. 


X 


Y 


J 


<rj 


H 


(Th 


K s 






(deg) 


(deg) 


(px) 


(px) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 



1 


14.79218 


-72.16765 


531.765 


155.835 


16.855 


0.013 


16.853 


0.012 


16.920 


0.023 


2 


I4.79437 


-72.16741 


548.261 


150.011 


16.011 


0.014 


15.968 


0.013 


15.979 


0.034 


3 


14.85005 


-72.16258 


965.954 


31.829 


17.483 


0.007 


17.098 


0.010 


17.097 


0.018 


4 


14.79519 


-72.16639 


554.363 


124.946 


17.957 


0.029 


18.110 


0.054 


18.113 


0.043 


5 


14.79583 


-72.16637 


559.162 


124.505 


15.665 


0.004 


15.664 


0.009 


15.731 


0.018 


6 


14.79584 


-72.16676 


559.283 


134.180 


18.141 


0.025 


17.840 


0.037 


17.714 


0.055 


7 


14.80519 


-72.16749 


629.400 


151.956 


18.724 


0.023 


18.307 


0.026 


18.074 


0.027 


8 


14.81704 


-72.16776 


718.270 


158.700 


19.061 


0.017 


18.626 


0.031 


18.561 


0.033 


9 


14.84279 


-72.16213 


911.515 


20.758 


17.483 


0.006 


17.117 


0.010 


17.097 


0.016 


10 


14.85175 


-72.16175 


978.702 


11.646 


19.181 


0.021 


19.364 


0.065 


19.252 


0.059 


11 


14.85197 


-72.16337 


980.320 


51.276 


18.321 


0.010 


17.905 


0.022 


17.907 


0.023 


12 


14.78242 


-72.16232 


458.576 


25.405 


17.969 


0.014 


17.820 


0.019 


17.878 


0.023 


13 


14.78620 


-72.16421 


486.962 


71.662 


16.760 


0.008 


16.605 


0.009 


16.645 


0.011 


14 


14.79341 


-72.16464 


541.046 


82.180 


18.041 


0.008 


17.567 


0.016 


17.488 


0.014 


15 


14.80191 


-72.16452 


604.786 


79.288 


18.263 


0.011 


17.814 


0.019 


17.679 


0.017 



of its newly discovered stars were determined. Approximately 
30% more stars were found in the region of NGC 346 from the 
second photometry run, but most of them have large photometric 
uncertainties. The final numbers of identified stars in each filter 
is 1 1 900 in J, 6 406 in H, and 5 837 in K s . We match these pho- 
tometric catalogs with a procedure developed in IDL to identify 
stars in common. After selecting the stars with the highest pho- 
tometric quality, i.e., stars with photometric uncertainties equal 
or smaller than 0.1 mag, we identified 2 783 stars in both the J 
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Fig. 5. Comparison of our photometry with that derived by 
2MASS for more than 20 common stars, located in all three ob- 
served fields for all three J, H, K s filters. Top: ISAAC versus 
(vs.) 2MASS photometry for the common stars. The linear fits 
that define the offset for each filter are shown with different lines 
and colors. The corresponding offsets applied to the absolute cal- 
ibration of our photometry are also given for every filter. Bottom: 
Dispersion in the 2MASS-ISAAC magnitudes for the common 
stars as a function of ISAAC magnitudes. The offsets for each 
filter are again defined by different lines. The standard deviation 
in the points derived from the corresponding mean values is the 
measurement of the uncertainties in our photometric calibration 
for each filter (see Sect. 13.21 ). 



and H filters, 3 067 stars in J and K s , 2 550 stars in H and K s , 
and 2 506 stars in all three filters collectively. 

3.1. Completeness and photometric accuracy 

Typical uncertainties in our photometry are shown in Fig. [2] The 
completeness of our photometry is evaluated by artificial star ex- 
periments, by adding artificial stars into the original frames with 
the DAOPHOT subroutine addstar. A few hundred artificial 
stars were added to each observed frame for every filter. This 
process was repeated several times, and the completeness of our 
photometry was derived by applying the photometric process to 
each new frame and comparing the final stellar catalogs with 
those of the original artificial stars based on the coordinates of 
the stars. The derived completeness of our photometry is shown 
in Fig. [3] for all three bands. Based on the completeness mea- 
surements and the photometric uncertainties, the limiting mag- 
nitudes of our photometry are 7 ^ 22.0 mag, H 20.5 mag, and 
K s si 20.5 mag. 

3.2. Photometric calibration 

We calibrate the magnitudes of the stars with near-IR absolute 
photometry of the same field from the Two Micron All Sky 
SurvejQ (2MASS). The 2MASS filter system agrees very well 
with that of ISAAC, and the 2MASS Point Source Catalog (All- 
Sky 2003) is therefore particularly suitable for the calibration of 
our photometry. The application ALADIN3 was used to identify 
all stars found with 2MASS in the same field-of-view as ours 
and to compare the magnitudes measured from both ISAAC and 
2MASS data sets. To calculate the magnitude offsets between 
the two photometries, we selected more than 20 bright stars with 
the highest quality photometric data from our photometric cata- 
log in each filter. The comparison of our photometry with that 
of 2MASS is shown in Fig. [5] The offsets per filter are cal- 
culated in terms of the median of the differences between the 
magnitudes of the stars, which are common to both ISAAC and 
2MASS samples. The apparent magnitude offsets are 1.812 mag 
in J, 0.678 mag in H, and 0.367 mag in K$. The uncertainties 
in this calibration can be quantified by the standard deviations 



6 http://irsa.ipac.caltech.edu 

7 http : // aladin . u- strasbg . fr/ aladin . gml 
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Fig. 6. CMDs of all stars found in all observed fields of the region NGC 346/N66 and its nearby control field of the SMC. CMDs in 
different combinations of the three bands are shown. Spec ifically, (a) J — H, H , (b) H — K s , K s , and (c) J - K s , K s . Two indicative 
isochrones from the Padova grid of evolutionary models (Gir ardi et al.| [2002), corresponding to ages ~ 4 Myr and ~ 4.5 Gyr are 
over-plotted to demonstrate the differences of the comprised stellar populations. 



derived from the dispersion in 2MASS-ISAAC magnitudes for 
the common stars (Fig. [5]- bottom). They are crj = 0.14 mag, 
cr H = 0.12 mag, and <jk s = 0.13 mag, respectively. 

3.3. The photometric catalog 

We derived celestial coordinates for all stars detected with high 
quality photometry (cr < 0.1 mag) from their pixel coordi- 
nates in the final jitter FITS images and used the applica- 
tions xy2sky and sky2xy available from the World Coordinate 
Systems (WCS) Tool$ (ver. 3.7.2). We transformed the (X,Y) 
coordinates of the stars into celestial coordinates according to 
the astrometric corrections provided in the FITS header of the 
corresponding jittered images by applying xy2sky. We then 
transformed the celestial coordinates of all stars with the use 
of sky2xy into a common (X,Y) pixel coordinate system with 
respect to the central field, where NGC 346/N 66 is mainly ob- 
served, by using the astrometric corrections provided in the FITS 
header of the K s image of this field. A sample of the final com- 
piled catalog of the stars found in all three filters is given in 
Table Q] This table is available in its entirety at the CDS. In 
Fig. HI the map of all sources detected in all three wavebands 
is shown. 



4. Observed stellar populations 

In our subsequent analysis, we consider only sources detected 
with the highest quality photometric parameters and with photo- 
metric uncertainties based on the PSF fitting of cr < 0. 1 mag in 
every waveband. The region NGC 346/N6 6 is known to h ost a 
mixture of stellar populations (e.g.. lGouliermis et al.ll2006h . in- 
cluding the evolved SMC field stars, the young main-sequence 

8 Available at http: //tdc-www. harvard, edu/wcstools/ 



(MS) and pre-main sequence (PMS) populations of the associa - 
tion NGC 346 and its vicinity (e.g.. iHennekemper et alj |2008). 
and the ~ 4.5 Gyr-old faint MS and bright RGB stars of the clus- 
ter BS 90 dRochau et al.ll2007l) . 

4.1. Color-magnitude diagrams 

The color-magnitude diagrams (CMDs) were compiled from 
our photometry for all stars detected in both the area of 
NGC 346/N66 (where BS 90 is also included) and the con- 
trol field, specifically the J-H versus (vs.) H, H - K s vs. 
K s , and J - K s vs. CMDs, shown in Fig. [6] The variety of 
stellar types in the observed region is illustrated by the super- 
imposed isochrone models. These evolutionary models, which 
are designed for the ESO Imaging Survey WFI UBVRIZ and 
SOFI JHK VEGAmag systems, were developed by Girardi et al. 
(2002]). The superimposed isochrones, corresponding to ages 
~ 4 Myr and ~ 4.5 Gyr, were selected to correspond to the 
metallicity of the SMC, namely Z = 0.004 dHaseret all 1 19981 
iBouret et al.ll2003l) . Both the bright MS of NGC 346 and the red 
giant branch (RGB) of the field and BS 90 are clearly defined in 
the CMDs of Fig. [6] these features being the most clearly distin- 
guished in the J-K S ,K S CMD. To determine the interstellar red- 
dening of the observed stellar populations, we assume a Galactic 
interstellar extinction law 

where Ry is roughly equal to 3.14 + 0.10 (ISchultz & Wiemeil 
fl975t) . The Galactic extinction in near-IR wavebands is given by 
A ., = 0.282 x A v A H = 0.175 x A v and A K = 0.112 x A v 
dRieke & Lebofskvlll985l) . Assuming these factors, we found 
that to fit the young isochrone to the observed bright MS pop- 
ulation in the CMDs of Fig. [6] an extinction of Ay ^ 0.3 mag 
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Fig. 7. / - K^, K s CMDs of stars found in two selected parts of the observed region NGC 346/N66, for (a) a circular area of radius 
~ 10 pc centered on NGC 346, (b) a circular area of radius ~ 8.5 pc centered on BS 90, and (c) the control field alone. The variations 
in the CMD from one area to another are highlighted by the superimposed isochrone models and typical positions of early-type MS 
stars and late-type giants (see Sect. 14.21 1. 



should be applied. On the other hand, the observed RGB stars 
do not appear to be affected by any significant reddening in the 
CMDs (Ay — 0.02 mag). This implies that the younger MS stars 
are still embedded within dust, while the older population of the 
field and BS 90 are not located in regions of high nebulosity, 
i.e., BS 90 is located in front of N66, and consequently their 
light does not suffer from extinction. 

4.2. Variations in the color-magnitude diagram 

The stellar systems NGC 346 and BS 90 represent two quite 
different types of clusters both from morphological and evolu- 
tionary point-of-views. Their evolutionary difference is naturally 
based on the stellar member populations of each system, which 
should define the prominent features in the corresponding CMDs 
of the systems. This becomes clearer when the CMDs of selected 
areas centered on these systems, comprising the most representa- 
tive stellar populations in the systems, are constructed. In Fig. [7] 
the J-K s , K s CMDs of (a) a circular area of radius ~ 35" (10 pc) 
centred on NGC 346, and (b) a circular area of radius ~ 0.5' 
(8.5 pc) centred on BS 90 are shown. These radii do not cor- 
respond to the sizes of the systems, but are selected to be very 
close to their centers to isolate the most prominent stellar pop- 
ulations in each system. The corresponding CMD of the control 
field is also shown for comparison (Fig. [7};). 

The variations in the CMD from one area to the next are indi- 
cated by the superimposed isochrone models in Fig. [7] The stel- 
lar population of the star-forming association NGC 346, as seen 
in the CMD of Fig. [TJ a), includes mostly young bright MS stars 
with a small contamination of the general field of the SMC, as 
this is defined by the CMD of the control field shown in Fig.|7fc). 
Typical positions of e arly-type stars are als o shown in Fig. |7Ja) 
based on the study of Hanson et al. (1997). On the other hand, 
the stellar content of the intermediate-age star cluster BS 90, 



seen in the CMD of Fig. |7jb), is characterized by the promi- 
nent RGB and few supergiants of age ~ 4.5 Gyr. The pos itions 
of late-type giants ar e also p lotted based on lPeletierl d!989l) and 
iBiimev & Merrifiejd] ([1998). Isochrone fitting on the CMDs of 
Fig-Hallows us to derive an indicative age for the stellar systems 
in the region. For NGC 346, we find an age ~ 3.9 Myr, consistent 
with that found by, e.g. jNota etall d2006l) . The CMD on BS 90 
is most accurately fitted by the isochrone of an age ~ 4.5 Gyr, in 
good a gree ment with the value p reviously found bv lSabbi et al.l 
d2007l) and lRochauetalJ d2007l) . It should be noted that BS 90 
is a populous cluster, of large spatial extent, and therefore the 
contamination of the stellar populations observed in the area of 
the association NGC 346 by those of BS 90 should not be neg- 
ligible. We discuss this in detail in Sect. 15.11 Finally, the CMD 
of the control field provides evidence of a variety of stellar pop- 
ulations, which are the results of different star formation events, 
of a wide range of ages, from young MS populations with age 
~ 30 Myr, an order of a magnitude older than those in NGC 346, 
to evolved giants of age > 2 Gyr. We note that near-R photome- 
try of early-type MS stars is quite insensitive to ages. Therefore, 
the ages given here are only indicative. 

The bright main sequence stars in NGC 346/N66 are known 
to be O-, B-, or early A-type stars of masses between 3 and 60 
M , while its low- and intermediate-mass stellar content is char- 
acter ized by a large sample o f pre-main sequence (PMS) stars 
(e.g jHennekemper et alj2008l) . The most massive of these stars, 
such as Herbig Ae/Be stars, are possibly also located in the red 
part of the main-sequence of Fig. EJa). However, our photome- 
try is not in general deep enough to detect significant numbers of 
low-mass PMS stars. The part of the CMD between th e MS and 
the R GB should also host classical Be stars (see e.g.. lBik et al.l 
2006). To perform a more accurate identification of these sources 
and to utilize the near-IR three-band color-color diagram of the 
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Fig. 8. Color-composite images of the main part of the region NGC 346/N66. Left: Image constructed with the combination of the 
jittered images in J (blue), H (green), and K s (red) from VLT/ISAAC imaging. Right: Image produced from the K s ISAAC image 
(blue) and Spitzer/IRAC 3.6 //m and 8 yum filters (green and red respectively). ISAAC images were constructed with DIMSUM 
(Deep Infrared Mosaicing Software), an external package of IRAF, especially designed to produce accurate sky subtracted images 
from dithered observations. 



detected sources in the young association NGC 346, we focus 
our subsequent analysis on the central area of NGC 346/N66. 

5. Stellar sources under formation 

There are no comparable studies to ours in the SMC, but there 
have been a few previous near-IR photometric investigations 
of resolved populations in star-forming regions of the low- 
metallicity environment of the Large Magellanic Cloud (LMC). 
A combination of ground-based near -IR data and space optical 
observations of the 30 Dor Nebula by Rubi o et alJ dl998h helped 
identify numerous stellar IR sources in or near the bright nebu- 
lar filaments west and northeast of R136, suggesting that a new 
stellar generation is being produced by the energetic activity of 
the massive central cluster in the remanent interstellar material 
around its periphery. Near-IR photometry of the second largest 
Hn region in the LMC, N11 B, also highlighted several bright 
IR sources dBarba et alJl2003h . Several of these sources have IR 
colors with YSO characteristics, and they are prime candidates 
to be intermediate-mass Herbig Ae/Be (HAeBe) stars. 

Deep near-IR imaging of the N 1 59/N 1 60 star-forming region 
in the LMC detected candidate HAeBe stars do wn to ~ 3 M , 
based on their near-IR colors and magnitudes (Nakaiimaet al. 
120051) . Two embedded HAeBe clusters were also discovered, one 
of them, and two neighboring Hn regions, providing hints of the 
beginning of sequential cluster formation in N159S. The spatial 
distributions of the HAeBe and OB clusters indicate large-scale 
sequential cluster formation over the entire observed region from 
N 160 to N159S. Near-IR photometry obtained to study the stel- 
lar content of the LMC star-forming region N4 is used to study 
the con nection between the interstellar medium and its stellar 
content (IContursi et al.l l2007). This analysis found several bright 
IR sources with characteristics of massive, early O-type stars. 
However, according to these authors, IR spectroscopy of these 
sources would confirm their very young and massive nature. 



IChen et ail d2009l) presented the most complete identifica- 
tion of YSO candidates in the LMC Hn complex N44. These au- 
thors combined mid- and far-IR Spitzer Space Telescope imag- 
ing with complementary ground-based imaging in UBVIJK to 
classify the YSOs into Types I, II, and III, according to their 
spectral energy distributions (SEDs). In their sample, ~ 65% 
of these objects were resolved into multiple components or ex- 
tended sources. The distribution of the YSO candidates com- 
pared with those of the underlying stellar population and inter- 
stellar gas illustrates a correlation between the current formation 
of O-type stars and previous formation of massive stars, provid- 
ing evidence of triggered star formation in N44. 

In connection to the aforementioned studies, our photometry, 
although not very deep, provides original near-IR colors of can- 
didate young sources in the vicinity of N66. However, from near- 
IR photometry alone it is impossible to accurately determine the 
nature of the most IR-bright sources, unless this photometry is 
part of a multi-band optical and IR investigation, as we discuss 
later. In the present study, we provide only a first set of accurate 
near-IR colors for all young objects in the region, and select the 
most interesting IR-bright objects from our photometric catalog, 
based on their near-IR colors and their near-IR excess inferred 
from the J - H, H - K s color-color diagram (C-CD). This selec- 
tion is described in the following section. 

5. 1 . Selection of young stellar sources 

To identify the sources that represent the most recent star for- 
mation in the region, we first consider the contribution of both 
the evolved MS and RGB stellar populations of the SMC field, 
the cluster BS 90, and the association NGC 346 to the com- 
plete stellar sample. As shown in Figs. [6] and [7] the MS is 
well defined with its brighter members belonging to the asso- 
ciation, but it is its red part and the region between the MS and 
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Fig. 9. J - K s vs. K s color-magnitude diagram (CMD) and the corresponding J - H vs. H - K s color-color diagram (C-CD) of all 
sources found in the main part of the star-forming region NGC 346/N66 with high photometric accuracy (cr < 0.1 mag) from our 
VLT/ISAAC photometry. Typical locati ons of Herbig Ae/Be stars (HAeBe; yellow shaded areas) and classi cal Be stars (red shade d 
areas) are shown based on the studies bv lLada & Adamsldl992l) . lDou ghertv et al. ( 1994). lMever et alj d 1997t) and lEiroa et al.| j2002). 
While T Tauri stars (TTS) are fainter than our detection limit in the CMD, typical positions of classical T Tauri stars (CTTS; blue 
shaded area), and the general loci of weak-line T Tauri stars (WTTS) are also shown in the C-CD for comple teness. In the CMD, 
the g reen "★" symbols represent the extinction-corrected positions of the massive YSOs from the samples of lHanson et alj d 19971) 
and lBik et al.1 (£006). The green and red lines in the C-CD denote the loci of MS and RGB stars, respectively. Reddening vectors 
for Ay = 5 mag and 10 mag are also drawn. These diagrams can be reliably utilized to distinguish the sources with infrared excess 
located to the right of the reddening vectors, away from the evolved MS and RGB stars. 



RGB that host stellar sou rces at their earlier stages of formation 
(see e.g.. iBik et alj|2006l) . Therefore, and to eliminate the sig- 
nificant contamination of this part of the CMD by the general 
field and BS 90, which is very close to the central part of N66, 
we constrain our analysis to the ISAAC FOV that covers the 
main body of NGC 346/N66, defined as the central field with 
-72.21 < Decl.(deg) < -72.15 in the map of Fig. H Color- 
composite images of this field, which covers the association 
NGC 346, the bar of N66, and several subclusters and Hn re- 
gions in its vicinity, are shown in Fig. [8] The dust emission, seen 
in the 8 fim band, indicates the centers of the most recent star 
formation and demonstrates the different information revealed 
by different wavelengths for the same region. 

The use of the near-IR C-CD is a reliable method for 
detecting sources, which are characterized by near-IR excess 
emission, from their positions in the diagram. In this dia- 
gram, the evolved populations are easily identified by compar- 
ison with models and considering interstellar extinction, PMS 
stars, namely classical and weak- line T Tauri (CTT, WTT) stars 
(e.g JAppenzeller & Mundtll989l). Herbig Ae/Be (HA eBe) stars 
(e.g. JPerez & Gradvll997U Waters & Waelkensll998h . and other 
YSOs (e.g.. lLada& Wilkinglll984HAndrelll994 are located in 
the reddest part. This is a signature of their IR excess due to their 
circumstellar dust in the form of cocoons or disks. Hence, when 
selecting candidate YSOs and PMS stars, we consider only the 
sources that are detected in the central ISAAC field in all three 
JHK S wavebands. The J - K s , K s CMD and J - H, H - K s 



C-CD of these sources are shown in Fig. [9] Considering that 
BS 90 to be a larg e cluster w ith a tidal radius bet ween r t 2! 15 
(ISabbi et alJl2007l) and 3: 15 dRochau et al.ll2007l) . its spatial ex- 
tent covers the central field of NGC 346/N66, and consequently 
affects the stellar content of the region. This is clearly demon- 
strated in the CMD of Fig. |9l where a prominent RGB can be 
seen. While these stars are bright in the near-IR, they do not ex- 
hibit any significant excess emission. Our classification to iden- 
tify young stellar sources based on their near-IR excess, there- 
fore allows us to effectively discard most of the evolved stel- 
lar contaminants from our sample of selected young sources in 
formation. However, a definite selection requires a sophisticated 
multi-band study on an individual source basis, as we discuss in 
Sect. 1331 

In the J - K s vs. K s CMD of Fig. [9] the areas oc- 
cupied by HAeBe (|Eiroa et alJ |2002) and classical Be stars 
dDoughertv et aljfl994l) are delineated by two ell ipses, and the 
positio ns of massive Galactic YSOs identified by lHanson et al.1 
d 19971) and Bik et al.l d2006l) are indicated by green stellar sym- 
bols. In the C-CD of Fig. |9j the typical positions of T Tauri 
(TTS), HAeBe, and classical Be stars are drawn to indicate the 
loci, where these sources in our sample should be expected in 
this diagram. We note, however, that typical TTS are fainter than 
the detection limit of our photometry and therefore are barely 
discernible in the CMD of Fig. |9l and only a few of them be- 
ing seen in the C-CD. We refer, however, to these PMS stars 
for reasons of completeness. The diagrams of Fig. [9] appear to 
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Fig. 10. Selection of the sources identified to have a near-IR excess based on their locations in the H-K s , J-H C-CD, and J-K s , K s 
CMD (see Fig . [9j» . Left: The fraction of the color indices (H-K S )/(J-H) as a function of the K s magnitudes of the detected sources. 
Right: The C-CD diagram of all sources selected as the most prominent young stellar candidates with near-IR excess emission. 



comprise a large variety of different stellar types at quite dif- 
ferent evolutionary stages. The regions in the CMD and C-CD, 
where classical Be stars are typicall y expected, also in clude MS 
populations. Searching for Be stars. iKeller et al J (1 19991) found in 
six fields centered on young clusters of the Magellanic Clouds 
(NGC 346 included), that the average fraction of Be to normal B 
stars is similar to that found in the Galaxy (~ 20%, see, e.g. JCoxl 
2000, p. 414). No connection between the Be star fraction and 
age or metallicity was found by these authors, and the classical 
Be stars detected in NGC 346 do not have any influence on the 
evolution of the region. According to their near-IR colors, red gi- 
ants and subgiants are expected, to occupy the blue edges of the 
typical locations of HAeBe stars. As a consequence, we should 
be careful to differentiate the sources that are most probably at 
their earlier stages of evolution, based on their near-IR excess, 
from the evolved red stars of the region. 

Based on the discussion above, we make a first tentative se- 
lection of the sources that most probably have near-IR excess us- 
ing a diagram of the color fraction (H-K s )/ (J-H) of all detected 
objects as a function of their brightness in K s . In this diagram, 
shown in Fig.[l0](left), one can see that most sources are located 
along a horizontal sequence of stars with (H-K S )/(J-H) 0.0, 
corresponding to the RGB and MS stars that exhibit no near-IR 
excess. Young stellar sources with a near-IR excess should be 
located away from the horizontal sequence. We first assume that 
all sources with (H - K s ) ~ 0.0 should not have any near-IR ex- 
cess, as they are located blueward of the MS in the C-CD of 
Fig. [9] We then apply a first-order selection of the remaining 
sources based on the criterion that an absolute offset from the 
horizontal sequence in the color fraction (H - K S )/(J - H) of 
about 0.3 is a reasonable limit to separate the stars with near-IR 
excess from those that show no excess. As a consequence, we 
select as sources with near-IR excess those that have color frac- 
tions (H - K S )/(J -H)> 0.3 or (H - K S )/(J - H) < -0.3 in the 
diagram of Fig. [TOl(left). 



However, an investigation of the positions of these near- 
IR bright sources in the CMD of Fig. [9] shows that there is 
significant contamination of these objects by classical Be and 
RGB stars, which should be eliminated. Therefore, we place yet 
tighter constraints on our sample by selecting only the sources 
that fall at the positions expected for HAeBe stars as shown in 
Fig. [9] (left) and redder. More precisely, we select the sources 
that fulfill the following criteria: (i) They are located to the blue 
part of the MS in the C-CD with (H - K s ) > 0.0. (ii) They 
have fractions of color indices (H - K S )/(J - H) > +0.3 or 
< -0.3. (Hi) They are located to the red part of a diagonal line 
that tangentially crosses the blue edge of the HAeBe area, spec- 
ified in the J - K s , K s CMD of Fig. [9] These sources repre- 
sent our final sample of objects selected as the most likely star- 
forming candidates, and the subject of our analysis here. The 
positions of these sources are indicated by red points in the dia- 
gram (H - K S )I( J - H) vs. K s of Fig. [10] Their positions in the 
C-CD are also shown in Fig. [10] (right panel). In this C-CD, it is 
indeed evident that our selected sample probably corresponds to 
objects with strong near-IR excess emission. 



5.2. The sample of candidate young stellar sources in 
NGC 346/N66 

The selection scheme returned 263 candidate young stellar 
sources in the main part of the region NGC 346/N66. The posi- 
tions of these sources are shown in a map of this area in Fig.fTTI 
as circular blue points. They are superimposed on a color- 
composite image constructed by combining our ISAAC image 
in the K s filter (green component) with ar chived image s obtained 
with the Infrared Array Camera (IRAC; iFazio et aUl2004l) on- 
board the Spitzer Space Telescope and in particular channel 4 
(8.0 yimi, red component), within the GTO science program 63 
(PI: J. R. Houck). These IRAC data were used to detec t candi - 
date YSOs in the general region of N66 bv [Simon et ail (120071) . 
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Fig. 11. Positions of stellar sources in our photometric catalog, with strong near-IR excess emission, based on their positions in 
the J - H, H - K s C-CD and J - K s , K s CMD (see Fig. |9). They are probably young stellar sources in the process of formation 
such as intermediate-mass PMS stars and high-mass YSOs. Their positions are overplotted with blue circle symbols on a color- 
co mposite image co nsisting of the K s ISAAC image (green) and SpitzerflRAC 8 /mi filters (red). Candidate YSOs from the catalog 
of lSimon et al.l d2007l) are overplotted with large red circular symbols. 



The objects identified by these authors as "definite" YSOs are 
also plotted in Fig. [TT] with large red circular symbols. 

In this map the positions of the detected candidate YSOs 
in NGC 346/N66 clearly correspond to dusty structures seen in 
the 8 /im emission as red filaments, these objects being located 
at peaks of mid-IR emission in these filaments. These peaks 
can also be seen - at lower spatial resolution - in the 24 fim 
band of Multiband Imaging Photom eter for SIRTF (MIPS; e.g. 
iHeim et alJll998t iRieke et al.ll2003l) onboard Spitzer, and they 
are also observed in the 2.14 nm H 2 line and the ISOCAM LW2 
band that covers 5-8 jum (IContursi et al.l 120001 : iRubio et"aT] 
2000). We note that the spatial distribution of our candidate 
young stellar sources follows the filamentary dust emission, be- 
ing concentrated along a few southern dusty arcs (at the middle 
and bottom of the map) and one northern dusty arm (at the top 
of the image). The latter, along with the south-eastern filament, 
is understood to be the product of a more recent triggered star 
forma tion event (Gouliermis et af] 120081) . iHenn ekemne r et al] 
(2008) discuss t he spatial distribution o f the low-mass PMS 
stars detected by Gouli ermis et al.l d2006l) . and of sources with 



excess Ha emission that are the most likely candidates to be 
intermediate-mass PMS stars (their Figs. 1 and 2). The spatial 
distribution of all these objects follows the same trend as our 
sources, and it is probable that both catalogs include several ob- 
jects in common. In addition, the positions of a significant num- 
ber of candidate young stellar sources found in our near-IR pho- 
tometry, shown in Fig. [TT] also seem to be quite clustered in 
several concentrations. We also note that almost all of these con- 
centrations coincide with subclusters of PMS stars observed with 
HST/ACS iden tified by applying two cluster analysis techniques 
dSchmeia et al.ll2009l) . 

We first compared a nominal search-box of 2" of the cata- 
log of our sources with that of the YSOs found with Spitzer by 
ISimon et al] d2007|) and PMS stars w ith Ha excess found with 
Hubble bv lHennekempereT al. (2008) in the same region. This 
comparison returns more than ~ 40% of the Ha excess stars and 
~ 60% of the YSOs also detected in the near-IR by ourselves. 
There are cases where more than one objects found with Hubble 
coincide with one of our sources, while in other cases a few of 
our sources are found as components of one of the YSOs found 
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with Spitzer within the specified search-box. This phenomenon 
is caused by differences in the resolving efficiency of the three 
instruments, demonstrating the importance of resolving multiple 
sources into their true components and classifying their true na- 
ture. The resolution achieved by VLT/ISAAC is at least 10 times 
higher than that of Spitzer/IRAC, allowing a reasonable iden- 
tification of any components in multiple YSOs. However, 15% 
of our sources coinciding with Ha-excess objects identified by 
iHennekemper et al.l d2008l) . are resolved by HST/ACS to be mul- 
tiple systems themselves. Therefore, the differences in resolution 
between the data sets obtained with HST, VLT, and Spitzer is 
a crucial issue in identifying true single-objects, or the compo- 
nents of multiple systems in our sample. As we discuss later, to 
perform a more detailed study in the near-IR, observations of the 
highest possible spatial resolution are essential. 

It is almost certain that our sample of candidate PMS stars 
and YSOs is incomplete, since for example TTS are not included 
because of our detection limit. Other young stellar sources such 
as HAeBe stars may also be missing due to the strong constraints 
of our selection criteria. In addition, our sample may be con- 
taminated by evolved stars, which in general do not exhibit any 
significant near-IR excess. While the positions of most of the 
sources in our catalog coincide with the dust filaments of N66, 
illustrating their youthfulness, there are candidate young sources 
in our sample that are located away from the dusty filaments of 
Fig.QT| These are usually assumed to be evolved stars, or even 
background galaxies, rather than stellar sources in formation. 
However, decontaminating our sample of old stars in the field 
and BS 90 located in the central region of NGC 346/N66 on a 
statistical basis for such a small sample, or based on their po- 
sitions away from the dust filaments, will possibly compromise 
the catalog of sources with true near-IR excess and it produce 
selection effects. 

A more sophisticated selection on a source-by-source ba- 
sis is certainly required to identify the most prominent objects 
that represent the most recent star formation in the region. This, 
however, can only be achieved through the excessive use of 
imaging in many different wavebands and/or spectroscopy (e.g., 
IChen et al]|2009l) . so that complete SEDs of individual sources 
can be constru cted and consequen tly their true nature accurately 
defined (e.g., Whit ney et al.l l2008). Observations of higher spa- 
tial resolution will allow the components of any unresolved 
bright sources in our sample to be recognized. Such a thorough 
analysis would certainly include the use of previous observa- 
tions of the region NGC 346/N66 from vari ous telescopes at 
different wavelengths including those of HS T jG ouliermis et al] 
2006; [Hennekemper et al.l2008l) and Spitzer (Bo latto et al.l2007t 



Simon et al. 2007), as well as new near-IR observations of 



higher spatial resolution and sensitivity. A preliminary study of 
the available data and the preparation of follow-up observations 
in the near-IR is currently being performed by ourselves. 

6. Conclusions 

We have presented a detailed near-IR photometric study with 
VLT/ISAAC of the star-forming region NGC 346/N66 in the 
SMC. We have used archival ISAAC imaging of the general 
area of N66, which includes the stellar association NGC 346, 
the intermediate-age cluster BS 90, and a southern empty control 
field of the SMC. We have performed photometry on images ob- 
tained in the filters J (1.25 yum), H (1.65 fim), and K s (2.16 fim) 
and derived a catalog of more than 2 500 stars detected in all 
three wavebands. The color-magnitude diagrams of these stars 
include a collection of different stellar populations, comprising 



the young MS stars of the association mixed with the RGB and 
old MS stars of BS 90 and the general field of the SMC, but also 
objects that are at their earlier stages of their formation. We se- 
lect the best PMS and YSO candidates in our sample on the basis 
of their positions in the near-IR color-color and color-magnitude 
diagrams. 

We focus the selection of these sources on the central field 
observed with ISAAC, which covers only the main part of the 
nebula N66 and the association NGC 346 to avoid any severe 
contamination of our sample with the evolved red stars of BS 90 
and the field. In this area, our photometry detected 1 174 stars in 
all three wavebands for which the near-IR CMD and C-CD are 
constructed. In these diagrams, the evolved stellar populations 
are mostly aligned along the sequences of RGB and MS stars as 
they are expected to be by the evolutionary models depending on 
the interstellar extinction, but certainly contaminate the sample 
of young stellar sources. The reason is that, while these sources, 
such as low-mass T Tauri stars, intermediate-mass Herbig Ae/Be 
stars, and massive YSOs, exhibit excess emission in the near-IR 
due to circumstellar dust, their positions in the CMD and C-CD 
do not necessarily cover the reddest part. 

Bearing this in mind, we make a tentative selection of PMS 
and YSO candidates with (H - K s ) > 0.0 (redder than MS and 
RGB) in the diagram of the color fraction (H - K S )/(J - H) 
as a function of the K s brightness. In this diagram, sources that 
have an excess in their near-IR colors are located away from the 
horizontal sequence of evolved stars. We select the sources that 
have an offset from the horizontal sequence of \A(H - K s /J - 
H)\ > 0.3 and we decontaminate the sample by selecting sources 
located in the area of the J - K s , K s CMD that HAeBe stars 
are expected to occupy and with redder colors. We consider the 
selected objects as the most probable candidates to be stars in 
formation. This selection delivers 263 candidate young stellar 
sources, which are located along the dusty filamentary structures 
of N66 seen in the 8 /mi emission from Spitzer and the 2.14 fim 
H2 line. 

Combining observations at several wavelengths to construct 
complete SEDs of individual sources is the most accurate means 
of establishing their true nature. Objects from our catalog of 
young stellar sources do indeed coincide with candidate YSOs 
detected with Spitzer and sources with excess emission in Ha in 
the region observed with Hubble. However, since a large amount 
of data per object is required for detailed SED studies, it is nec- 
essary to enhance the available data sets with new data, prefer- 
ably obtained with cameras of higher resolving power, so that 
multiple objects can be resolved in their components. 
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